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Abstract. We classify multidimensionally consistent maps given by (formal or conver¬ 
gent) series of the following kind: 

00 

~ ^ij ”6 Xj ^ik' ' 

m=2 

where are homogeneous pol 5 momials of degree m of their respective arguments. The 
result of our classification is that the only non-trivial multidimensionally consistent map 
in this class is given by the well known symmetric discrete Darboux system 

rj, ^ij “6 ^ik^jk 

^k^ii ~ / , 

\/^ “ ~ ^jk 

1. Introduction 

The goal of this paper is to contribute to the problem of classification of discrete three- 
dimensional (3D) integrable systems. 

The notion of integrability we adhere to is the multidimensional consistency. It originates 
in the development of integrable structures of Discrete Differential Geometry, see Q for 
details and for historical remarks. Multidimensional consistency of a given discrete d- 
dimensional system means that it can be imposed in a consistent way on all d-dimensional 
sublattices of a {d + l)-dimensionl lattice. In retrospect, this property can be under¬ 
stood as a discrete analog of the existence of commuting hierarchies of integrable sys¬ 
tems. It has been shown that this property guarantees the existence of such integrability 
attributes as the discrete zero curvature representation and Backlund transformations en¬ 
joying Bianchi-type permutability. Moreover, these attributes are encoded in the system 
itself and can be found in a straightforward and algorithmic way. Thus, multidimensional 
consistency can be accepted as a rather general and powerful definition of integrability of 
discrete systems. 

In Ref. [[l|, a classification of 2D systems on quad-graphs with scalar (complex) fields 
at vertices, integrable in the sense of 3D consistency, was performed. As a result of this 
classification, a rather short but exhaustive list of integrable systems of this kind was 
produced, known nowadays as the ABS list. Moreover, one class of multidimensionally 
consistent 3D systems has been classified, namely those of the octahedron type 0. It was 
demonstrated that the relevant combinatorial structure in this case is a multidimensional 
root lattice of type A. The most striking feature is that the number of integrable systems 
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drops dramatically with increasing dimension: only half a dozen of discrete 3D systems 
with the property of 4D consistency are known, almost all of them are of a geometric 
origin. 

In the present paper, we classify 3D maps within certain general Ansatz, which are 4D 
consistent. The main result of our classification is that, within the considered class, the 
only non-trivial 4D consistent map is the symmetric discrete Darboux system (equation 
^ below). Note that this system appeared in the study of multi-dimensional circular 
nets 1^, see also 0, and enjoys several distinct geometric interpretations, for which one 
can consult Q and Q. 


2. The 4D consistency property 

Consider a 3D map <l>(x,i/,z) = (x^y,z). A possible definition of integrability of such 
maps is their 4D consistency Q, which we describe now. 

Combinatorially assign the quantities x = X 23 / y = z = Xi 2 to the three faces of 
a 3D cube parallel to the coordinate planes 23, 13, 12, respectively. Let stand for the 
unit shift in the k-th coordinate direction. We assign the quantities Xij = T^xij to the three 
opposite faces. Therefore we view the map O as 

‘hi23 (■T12/^13/'T23) = (’f3'Tl2/^2^13/7^1^23)/ (1) 

where the indices 123 of O represent the 3D space where the map acts. More generally, 
we assume that in any 3D sublattice of spaimed by the directions i,], k {i < i <k) 
there acts a map 




Figure 1. A map on an elementary 3D cube with fields assigned to 2-faces. 

Now consider the initial value problem with data Xij, i,] = 1,2,3,4, i < j, prescribed at 
six squares adjacent to one common vertex of the 4D cube. Then the application of a 3D 
map like Q to the four 3D cubes adjacent to this vertex allows us to determine all T]^Xij: 

‘hl23 (^12/'T13/^23) = {T3Xi2,T2Xiq„TiX23), 

^124 (t 12 /Ti 4 ,X 24 ) = (T4X12, T2X14, 74X24)/ 

‘hl34 (xi3,Xi4, X34) = (74X13,73X14,74X34), 

‘^234 (^23/'T24/T34) = (74X23,73X24,72X34). 
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At the second stage, the map is applied to the other four 3D cubes of the 4D cube: 

‘hl23 (14X12,14X13, T4X23) = (14(13^12)/14(12 ^13)/l4('ri^23))/ 

‘hl24 (l3'^12/13^14/13^24) = (1"3(r4Xl2)/13(12 X14), 73(71X24)), 

‘hl34 (T2X13, T2X14, 72X34) = (72(74X13), 72(73 X14), 72(71X34)), 

‘^234 (7 iX 23, 7 iX 24, 71X34) = (73(74X23), 7i(73 X24), 71(72X34)). 


Now, 4D consistency of means fhaf the following six equations are identically satis¬ 
fied for arbitrary initial data: 


TiiTkXij) = Tk{TeXij), i < ;, {z,;,;c,£} = {1,2,3,4}. 


( 2 ) 


In this setting, the two most celebrated 4D consistent maps are the symmetric discrete 
Darboux system, 

Xij Xi](Xj]( 

TkXij = 


1-4 


1-4 


(3) 


(where one assumes that x^ = xy/), and the star-triangle map. 


Xij 


^ij^jk 7 ^jk^ki 7 


(4) 


(where one assumes that x^y = —Xji). 


3. Statement oe the classieication problem 
We consider maps Q given by (formal or convergenf) series of the following kind: 


7 fcX;y X;y-|- (X/y, X;^, Xyj;), 

m=2 


(5) 


where are homogeneous polynomials of degree m of their respective arguments. 
Each time when the double indices are associated with a 2D coordinate plane (i.e., in 
notations like X/y, A^) we assume that the two indices are permutable, that is, X;y = Xy,, 

= Ay..^, etc. Note that the Darboux system (|3|) belongs to this class, while the star- 


ii;k' 


Xj 


triangle map ffl does not for two reasons: first, the ^ew-symmetry assumption x^y = — 
is enforced, and second, if is not close to identity. 

We will classify 4D consisfent maps of type We perform the classification modulo 
the group of admissible transformafions generated by the point (coordinate-wise) trans¬ 
formations 


Xij ^ Xij + Y^b 


m=2 


(m) m 
ij ^ij> 


given by (formal or convergenf) series, and by the scaling transformations 


Xij I }■ CijXij. 


( 6 ) 


(7) 


Here and c^y are arbitrary constants. Indeed, transformations from the group gener¬ 


ated by 0 and Q leave the form of <l>;yj- invarianf. 
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( 2 ) 

Our first aim is to determine quadratic polynomials which can appear as leading 

terms in 4D consistent maps satisfying (j^. All lemmas in the present section are devoted 

to such a characterization. 

(2) J I 

Lemma 1. Polynomials G {1,2,3,4}) o/4D consistent maps {5^ satisfy the following 

system of six nonlinear partial differential equations (in which {i,i,k,i} = (1,2,3,4}); 


a 4 (2) 

_^4(2) , 4(2) , A2) ^ 

dxii ax,.;, dxjj, AA 


aA; 


( 2 ) 




Wt (2) A2) Al) 

ax,-; ax,.^ axy^ AA- 




( 8 ) 


Here and below, if the arguments of are not written explicitly, they are supposed to be 

{xq, Xi](, Xjf). 

Proof. We compare the terms of degree 3 in Tg{T]^Xij) = T]^{TiXij). We have: 

2Cij + A--,j^(^Xij,Xi](,Xj]f) + A-_^(^Xij, Xi]^, Xjf) + . . ., 

T^xq = xq + A\J^^{xq, xu, Xji) + A^^^^{xq, xu, Xji) + ..., 

so that, taking into account only terms relevant for the expansion up to degree 3, 

Ti(T,xq) = xq + A^^ + A^, + ... 

+ ^ifk + ^ife + ■■■>Xik + a|^?^ + ...,Xjk + A^l^^ + ... j 

A ^Ipk + • • • / + • • • / + • • •) + • • • 

= Xq + aJ^?^ + + a|^?^ + a[^^^ 


+ 


A2) 

aX;; aA ^Xik 




34 ( 2 ) 

ijAM) ^^iyk A2) 


and similarly, 

T'i{T'kXq) = 2 Cq + A^-^^ + a |^..]^ + ... 


+ {^ij + + • • • / + aJ|\ + ..., xy^ + Aj|\ + ... 


A i^h + • • • /+ • • • / + • • •) + • • • 

= + ^§k + 4 a 


+ 


34 ( 2 ) 34 ^^^ 34 ^^^ 

_^4(2) , hA A2) .(2) 

dxq hA^ dxii 3x.^ ^y£;fc + ---- 


Now equations (|^ follow directly by comparing the terms of degree 3 in Tg{T](Xq) and 


TiiTj^xq). 


□ 
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Next, we analyze the system for polynomials We use the following notation 
for their coefficients: 

Xik, Xjk) = Xij. kXikXjk + ^f.k^ijXik + + Ay. kX^j + }i\p^k^jk + (9) 

Lemma 2. For polynomials {i,j,k G {1,2,3,4}) of the form (91 satisfying system (8>, we 
have 

Xij;£^ij;]c 0 / ( 10 ) 

for all {i,j,k,i} = (1,2,3,4}. 

Proof. We substitute 0 into the set of differential equations 0 thus obtaining a system of 
cubic polynomial equations which must be identically satisfied. The monomial xyXigXji 
can only come from the first term on the left-hand side, and the vanishing of its coefficient 
results in condition ( 10 ). □ 

We now deliberately choose to consider only the main two branches of conditions ( [10] >, 
namely 

(I) Xij.^k 7 ^ 0 arid ^ij;k = 0 lor all i,j,k G ( 1 , 2 ,3,4}. 

(II) O-ij-^k = 0 and Ay.jt 7 ^ 0 for all i,j,k G (1,2,3,4}. 

Accordingly we have the claims contained in Lemmas]^ and 

Lemma 3. In case (I), polynomials {i,j,k G (1,2,3,4}) of the form (91 satisfy system (8Hf 
and only if 

^ik^jk 

Xy-k 


efi = 0- 




(11) 




for all i,j,k G (1,2,3,4}, where Cy are arbitrary constants. Modulo gauge transformations, we 
can assume 


4 (2) _ Y Y 

■^y-k — -^ik-^jk 


for all i,j,k G (1,2,3,4} 

Proof. Assume ay-k 7 ^ 0 and \iyk = 0 for all i,j,k G (1,2,3,4} and substifufe 0 into 
the set of differential equations 0. Collecting coefficients of monomials XigXjgXik, we find 

- kilo = 0/ which gives 

kfk = kiln ( 12 ) 

for all {i,i,k,i} = (1,2,3,4}. Similarly, collecting coefficients of monomials XigXjgXkg, we 
find Xy.g{f>f^.k + ^ft,k) = which gives 

= -^f,k (13) 

for all {i,i,k,i.} = (1,2,3,4}. Combining conditions (12) and (j^ gives flf^.k = 0 for 

all i,],k G (1,2,3,4}. Collecting coefficients of monomials xj^Xjk, we find y\k-e‘^iiA ~ 9, 

which gives the remaining conditions p^^-^k = 0 for all i,j,k G {1,2,3,4}. We are now leff 
with 


(2) 

■dy;)c Xy-kXikXjk- 


( 14 ) 
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Substituting this into (|8|) we find the following conditions for the coefficients 


^ik;£^ii;k ^ji}k^ij-,£ 




^ik}£^ii;k ^i£-,k^ii;£- 


( 15 ) 


To solve conditions (15), we re-write them as 

^tj;£ _ <^ik;£ _ ^jk;i 
^ij-fk ^i£;k ^i£;k 

The second equality here implies: 

2 

^ik;£^i£;k ^jk;£^j£;k • ^k£' 

since this quantity does not depend on i, j. As a consequence, we find: 

r2 /I 

^2 ^ij;k^ik}i ' ^ii;k^ik;j ^ik^jk 

^ij-,k ~ rmr” ~ ' 


^jk} i^ik;i 




which coincides with = ‘^ik<^ik/ Cjj from (111. It is immediately verified that this neces¬ 
sary condition is also sufficient for ( [15] |. The gauge transformation (j^ brings a map with 
the quadratic terms ([T4]| with = Cij^Cjj^/Cij to a similar map with = 1. □ 


( 2 ) 

Lemma 4. In case (II) polynomials of the form (9) satisfy system (|8|) if and only if 

Pii;k ~ rij-,k ~ 


SO that 


for all i,j,k G {1,2,3,4}. 


4 (2) _ y2 


Proof. Assume = 0 and f 0 for all i,j,k G (1,2,3,4} and substitute (|^ into 
the set of differential equations (8|. Collecting coefficients of monomials we find 


^ij-jflpk ~ which gives = 0 for all i,j,k G (1,2,3,4}. Next, collecting coefficients 

of monomials xjgXp, we find ^ij-,kh^ipi ~ ’^hich gives the remaining conditions = 0 

for all i,j,k G (1,2,3,4}. We are left with = Xij.^x^j, which satisfy (8| without any 
further conditions on the coefficients D 


j(0 - 


4. Classification for case (I) 

The main result of our paper is contained in the following Theorem. 

Theorem 1. Any 4D consistent system of map {i,j,k G (1,2,3,4}) of type ([^ with 


^ii;k ~ ■^ikXjk 


is equivalent to the symmetric Darboux system (|^ modulo admissible transformations of type ([^. 
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Proof. For any rn > 2, the terms of degree m + 2 in the 4D consistency condition (j^ yield 
the following recursive system of equations for the polynomials 

^4('«+l) 

^ik^jk-j +^ik^it,i 3 ^.. + ^a^ki 3 ^.^ +XjeXM (16) 

^i£^ie;k ^ji^it.k ^ik^jk ^jk^ke ' 

(s) 

where F is a polynomial function of the with s < m. This is a system of six 

partial differential equations for twelve unknown polynomials , which can be re¬ 

formulated as a system of linear algebraic equations for the coefficients of these polyno¬ 
mials. We know that ( [T^ admits at least one particular solution (corresponding to the 
symmetric Darboux map), and we are going to study the general solution. It is obtained 
by adding an arbitrary solution of the corresponding homogeneous system: 

.(m+1) , .(m+l) , ij-rk , ^^ij-.k , ^^ij-.k ... _ 

^ik^jt,£ +^ik^ik-,£ +^ie^j£ 3 ^.. + Xi£Xk£ + Xj£Xke ^ (17) 

^i^-^i£;k ^i^^i£-,k ^ik^jk ^ik^ki ^ ^ — 0 - 

Lemma 5. The general solution of the homogeneous system p7| ) is given by 

= Vjk (jnhijxfr'^ - , (18) 

where by are arbitrary constants. 

It turns out that this freedom is exactly compensated by means of admissible gauge 
transformations: 

Lemma 6 . Given maps {i,j,k G {1,2,3,4}) of type ([^, a change of variables 

Xy HA Xy + byXy 

f S1 (7Tl I X1 

with m > 2 does not change the polynomials with s < m, while the polynomials A^,^ ’ 
get shifted by 

XikXjk [-mbyxT-^ + bi^xf^-^ + bjkxjf^^ . 

Theorem]^ is an immediate consequence of Lemmas]^ and□ 

It remains to prove Lemmas and The second one is proved by a straightforward 
computation. 

Proof {ofLemma^. We consider in ( [T7| terms of a low multidegree with respect to x^, Xj£, 
Xy£, belonging to mree distinct classes: 

(a) Terms of multidegree (1,1,0) with respect to x^, Xj£, X]^£, and of degree m — 1 with 
respect to xy, Xj^^, Xj^. 
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(b) Terms of multidegree (1,0,1) with respect to Xj£, and of degree m — 1 with 
respect to Xy, Xjj^, Xj^. 

(c) Terms of multidegree (0,1,1) with respect to xig, Xj£, xj^i, and of degree m — 1 with 
respect to Xjj^, Xj^. 

Terms belonging to the class (a) can come from the following part of ( [T7| >: 

Ox-- ^ik^]k~ 




dXij 


which, upon cancelation of the common factor x^Xj^, leads to 


OA 


m+l) 


ii;k 


dXi 


Pixfk + pix’ll + PsXikXjkX”- ^ 


(19) 


where 


Pi is the coefficient by Xj£X^ in Aj^. 


m+l) 


(m+l) 


P 2 is the coefficient by x^x^ in A-^.^ , 

P 3 is the coefficient by Xi£Xj£X^~^ in /dXy. 

Terms belonging to the class (b) can come from the following part of ([T7|>: 


^jk'^'ik-fi ’ Xi£Xk£- 


(m+1) 


aA 


m+l) 


ii;k 


^ .(m+l) 

.(m+l) ^^ii;£ 

XU^jf-k XikXki XjkXkg- 


dXik dXil 

which, upon cancelation of the common factor XnXki, leads to 


dXj£ 


^ . (m+l) 


dx 


ik 


-p^XjkxJl ^ + Psxfk + PeXikx"^ ^ + P7Xjkx”l ^ 


where 


P 4 is the coefficient by Xi£Xk£X^~^ in A 
P 5 is the coefficient by in 
P6 is the coefficient by X;^x”)“^ in BAj^/^Vax/^, 


(m+l) 


t] - ii;i 


• p 7 is the coefficient by x i£X^ ^ /dxj£. 

Finally, terms belonging to the class (c) can come from the following part of ([T7|>: 


7) A ^ 

.(m+l) , ii;k 

^ik^jk-J Xj£Xk£ 


(m+l) 


A 


. (m+l) 

— Xj£/i^^.^ — XikXki 


ax,Y 


dxi 


i£ 


^i£ 


which, upon cancelation of the common factor Xj£Xk£, leads to 


dXjk 


-PsXikXf^ ^ + P9X'ik + PwXikX^ ^ApilXjkXfj ^ 


where 


ps is the coefficient by Xj£Xk£xJ^ ^ in 


( 20 ) 


( 21 ) 
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• p 9 is the coefficient by in 

• pio is the coefficient by Xj£X^~^ in /dx^, 

• pii is the coefficient by XjiX^~^ in /dxj£. 


Cross-differentiation of expressions in 


(19), (20), leads to: 


0xij0xi]( 


3xij3xj]( 


dXi^dXj^ 


mpix"^ ^ + P3Xj]cX^- ^ 

{m - l)p6XikX^~^ +{m- l)p7Xjkx’ll~^, 

mpixff'^ + p3Xikx'l-~^ 

{m - l)pioXikX^-~^ + (m - l)pnXjk^^~^, 

^ ^m—l I ^m—1 I ^ ^m—1 

-P4% +nip5X-j^ +P7Xij 

^ ^m—1 I ^m—1 I ^ ^m—1 

~PsXjk +7npgX-k + PwXjj . 


Thus, we find the conditions 


Pl = P2 = P6 = Pll = 0, 

P3= {m- l)p7 = (m - l)pio, 

P4 = -mp9, p8 = -mp5. 

Taking into account the above conditions we can write our preliminary results as follows: 




''PA 

= (m-l)pioX/fcXyfcX'y ^ 

(22) 

dxij 


= mpgXjkX’^-'^ + psxfk + pwXjkX^-^, 

(23) 

Oxik 

i];k 

= mp7,XikX^k^ + pgX^k + PwXikX”-~^, 

(24) 

Oxjk 


where 

• (m — l)pio is the coefficient by Xi£Xj£X^~^ in /dxij, 

• pio is the coefficient by Xi£X^~^ in /dxji, and at the same time the coeffi¬ 
cient by Xj£x'll~^ in /dxii, 

• p 5 is the coefficient by in 

• —mp 5 is the coefficient by xy^Xjt^x^”^ in 

• p 9 is the coefficient by xj-^xj^ in A^J!k^\ 
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• — mp 9 is the coefficient by ^ in 

Now from equations (22)-(24) we find: 


(m+l) 


A 


(m+l) _ 


iy,k = V5Vj{ + lp9Xik^ik + PwXikXjkX';- 


m—1 


(25) 


where 


m—1 

ij 


in 


• pio is the coefficient by 

• —mp 5 is the coefficient by in 

• —mpg is the coefficient by XigX]((x'^~^ in 

From the characterization of pio it follows that this coefficient depends on i,] only, being 
independent of k,L Thus, we can set pio = mhij. Then the characterizations of —mp^ and 
—mpg yield that ps = —bjj^, pg = This allows us to finally re-write ( [2^ in the form 
(181, which proves Lemma □ 

5. Classification for case (II) 

In the case (II) we only have trivial 4D consistent maps. 

Theorem 2. Any 4D consistent system of map £ {1/2,3,4}) of type ([^ with 

4 (2) _ pt , y2 

^ij;k ~ ^^i];kXij, 


^ii;k 0/ equivalent to 

T'k^ij ~ fii',k{^ij)/ 

where fij;]^ are univariate functions such that 

fir,k°fiiA = fir,t° fiyk 

for all i,j,k,i G (1,2,3,4}. 


(26) 

(27) 

(28) 


Proof. We prove by induction that in the expression (5) for Tj^xy each polynomial of 

degree m is in fact a monomial xf- depending on Xy only By assumption, this is true 
for m = 2. We fix m > 3, assume the said property is fulfilled for all degrees less than m 
and derive the same property for degree m. One computes that the terms of degree m + l 
in Ti{T],xy) - T^{TiXy) are given by 


2Ay.kXyA^I + 


dA) 


dx 


hA , 2 , ijA . 2 , 

\y. iXy + C+k + 


dAf 


n 


ik 


dXjk 


h'A , 2 

^jk+^jk 


fm) 


-2Xy,,XyA^:l 


dXy 'J 


dA) 


(m) 




dx 


i£ 


a/t'”] 

-V 2 UA ^ 2 , m 

^it,k+i - + V^y 


m+l 


ji 


(29) 


where the term px^+^ comes from A^{/,, A^{/. with s < m. Now it is easy to show that 


,(®) 4(s) 


A 


[m) 


ij; k' ij; £ 


•y.j. does not depend on x/j-, Xjj^. Indeed, if the degree of A^J!^ with respect to would 

be d G [1, m] (with a non-vanishing leading coefficient), then the degree of the expression 
( [2^ with respect to Xf^ would be d -|- 1 (also with a non-vanishing leading coefficient), so 
this expression could not vanish identically 
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We conclude that any 4D consistent system of maps G {1,2,3,4}) of type 

g with ( [2^ is equivalent to ( [27[ > for some functions The 4D consisfency condition is 
m nothing but that the requirement that the functions fy-]^ commute as in (281. □ 


6. Conclusions 

The present paper has been devoted to the classification of 3D maps of fype ^ which 
are 4D consistenf. The resulfs of our classification problem are contained in Theorems 
and 1^ The most interesting finding is that the symmetric discrete Darboux system ^ is 
the only 4D consistent map belonging to our Ansatz. 

Future research will be devoted to the classification of 4D consistenf 3D maps which 
are nof a perturbation of the identity (thus including the star-triangle map (|^), as well as 
to extension into higher dimensions. Let us mention that no 4D map which would be 5D 
consistent is presently known. 
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